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Abstract 

Detailed knowledge about equilibration processes is of interest for various fields of physics, in¬ 
cluding heavy ion collision experiments and quantum quenched condensed matter systems. We 
study the approach to equilibrium at late times within two types of strongly coupled thermal 
systems in 3 -|- 1 dimensions: systems in the presence of (i) a non-zero charge density, or (ii) a 
magnetic field at vanishing charge density. Utilizing the gauge/gravity correspondence, we map 
the aforementioned problem to the computation of quasinormal frequencies around two particular 
classes of black branes within the Einstein-Maxwell theory. We compute (i) the tensor and vec¬ 
tor quasinormal modes of Reissner-Nordstrom black branes and (ii) the scalar, as well as tensor 
quasinormal modes of magnetic black branes. Some of these quasinormal modes correspond to the 
late-time relaxation of the above systems after starting with initial pressure anisotropy. We provide 
benchmarks which need to be matched at late-times by all holographic thermalization codes with 
the appropriate symmetries. 
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I. INTRODUCTION 

Systems far from equilibrium have been largely inaccessible for a long time despite being 
abundant in nature, arising in both high energy particle and condensed matter physics. 
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FIG. 1: Schematic plot of QNM sets. Two distinct sets of QNMs of h^y fluctuations around 
Reissner-Nordstrom black branes are shown. Set 1 (black dots) is the set which was found already in 
the case of a black brane with vanishing charge [I] . Set 2 (red triangles) is a set of QNMs wandering 
up the imaginary frequency axis as the black brane charge is pushed towards its extremal value. 
Thus, at large charge the purely imaginary QNMs of set 2 will dominate the late-time behavior of 
the dual field theory plasma. 

Experimental examples of such systems are the quark gluon plasma state generated in heavy 
ion collisions ElE], and quantum quenches of condensed matter systems nisi, for a review 
and references see e.g. [6]. The past few years have seen an increased interest in systems 
far from equilibrium. This renewed interest was spawned by the development of methods 
allowing the study of such systems at strong coupling [MQ] utilizing the gauge/gravity 
correspondence m- This correspondence translates the problem of a strongly coupled 
system far from equilibrium into the equivalent gravitational problem of Ending the metric of 
a time dependent geometry. This task, in general, involves solving the Einstein equations as 
a system of partial difierential equations with a particular set of initial conditions. The main 
question to be answered within this framework is: what exactly is the full time dependent 
process of equilibration? 

At late times, the question “How does the system approach equilibrium?” is identical to 
the question “How does the system respond to small perturbations around equilibrium?”. 
Depending on the symmetries and details of the gravitational setup, and on the symmetries of 
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the initial conditions, the aforementioned far from equilibrium systems generically asymptote 
to a particular equilibrium state at late times. In many cases this equilibrium state is 
going to be some kind of black brane geometry uni [12]. In other words, at late times the 
fully dynamic nonlinear problem turns into a linear problem of perturbations around black 
branes. Therefore, one can study the late-time behavior of far from equilibrium systems by 
considering perturbations around particular black branes, i.e. by computing the quasinormal 
modes of those black branes. Hence, comparing far from equilibrium results at late times to 
quasinormal mode results serves as a powerful consistency check of these highly non-trivial 
calculations. 

Quasinormal modes (QNMs) are the “eigenmodes” of a black brane geometry. Classical 
black brane horizons are capable of absorbing energy, but can not emit energy. In this 
sense perturbations around a black brane can lose energy into the horizon, and are thereby 
a non-conservative system described by a non-self-adjoint operator. Normal modes of a 
self-adjoint operator have real-valued frequencies, whereas quasinormal modes associated 
with a non-self-adjoint operator are complex-valued. Their imaginary part encodes the 
dissipation of a mode oscillating with the frequency given by their real part. Using the 
gauge/gravity correspondence the quasinormal modes have been identihed with the poles 
of retarded propagators of the operator dual to the relevant gravity perturbation [T3|. For 
example, the shear perturbation h^y of the metric is dual to the xj/-component of the energy 
momentum tensor T^y within the dual held theory. Therefore, the location of the poles of 
the retarded propagator {T^y^Xi) Txy{x2)) are identical to the quasinormal mode frequencies 
of the dual gravity perturbation h^y 

Numerous results are known for quasinormal modes of black holes, see for example HI, 
but there are some important gaps in the literature. Hence, in the present paper, we 
study quasinormal modes of metric (and gauge held) perturbations around two distinct 
black brane systems: the Reissner-Nordstrom black branes and magnetic black branes^. By 
magnetic black branes we refer to the solutions of Einstein-Maxwell-Chern-Simons theory 
in the presence of a constant magnetic held found and analyzed in [T5lfT8] . These magnetic 
black branes are dual to a (3 -|- l)-dimensional held theory at nonzero temperature and with 

^ Note that we distinguish between black holes with a compact horizon and black branes which have a 
planar horizon. In the present paper, however, we are exclusively concerned with black branes. 
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a nonzero magnetic field. For example, the analysis in the present paper applies to TV = 4 
SYM in a magnetic held for a f/(l)-subgroup of the R-charge. In [T3] - [T^ such setups were 
studied with an eye on applications to condensed matter systems in magnetic helds on one 
hand, and applications to heavy ion collisions on the other ng. Quasinormal modes of this 
system have not been computed directly before. 

In contrast to that, the Reissner-Nordstrom (RN) black brane solution in Anti-de Sitter 
space (AdS), corresponding to A/" = 4 SYM with a chemical potential and charge density 
for the R-charge, are well studied. However, while many QNMs have been computed for 
the AdSi black holes and black branes mi EDI EH, the AdS^ Reissner-Nordstrom black 
brane QNMs are not as readily available. Stability of AdS^ RN black holes against metric 
and electromagnetic perturbations was shown in [22] and [23], establishing linearized mas¬ 
ter equations for perturbations with vanishing spatial momentum. RN black branes are 
invariant under SO{3) spatial rotations. Hence, their quasinormal modes (at vanishing mo¬ 
mentum) can be classihed into scalar, vector, and tensor modes because the corresponding 
perturbations transform as a scalar, vector, or tensor under SO {3) spatial rotations. Some 
of the AdS^ RN black hole and black brane QNMs have been calculated before mEo], in 
particular the RN black hole scalar QNMs have been studied in [23] using the formalism 
developed in [23]) while the RN black brane vector QNMs were studied in [25], and the 
asymptotic QNMs (those which are inhnitely damped) were discussed in [26]. While [25] 
studies vector QNMs exclusively, [2Z] considered a 5-dimensional (single charge) STU black 
hole and hnds an analytic solution for the tensor fluctuation in the hydrodynamic limit, 
where frequencies and momenta of the perturbation are small compared to the temperature 
of the black brane. Similarly, in [28] such an analytic solution is found for the tensor fluctua¬ 
tion of the single R-charged RN AdS^ black brane, and in addition the lowest hydrodynamic 
vector QNM is computed. None of the tensor QNMs are within reach of the hydrodynamic 
approximation in this case. Furthermore, within the hydrodynamic limit, also the lowest 
lying scalar [29] , vector and tensor BOI QNMs have been studied. A Chern-Simons term has 
been added to the Einstein-Maxwell action in 011 and influences the fluctuation equations 
and hence the transport effects, in particular exhibiting chiral transport in the vector fluc¬ 
tuations. However, the non-hydrodynamic tensor and scalar QNMs of the RN AdS^ black 
brane have not been computed, as far as we know. There seems to be no systematic study 
of the metric tensor fluctuation modes at sizable momentum i.e. beyond the hydrodynamic 
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limit. Neither are we aware of a systematic QNM study of the Reissner-Nordstrom black 
branes over the whole range of charge densities, which runs from zero to the extremal value, 
Qexi at which the black brane has vanishing temperature. In this work we provide a study 
of the metric tensor fluctuation quasinormal modes at sizable momentum and for a large 
range of charge densities from zero to over 90 percent of the extremal value. 

Two Mathematica [32] notebooks, see [33] for download and details, accompany this pub¬ 
lication, allowing the reader to compute quasinormal modes -with any desired accuracy- for 
any given charge density and momentum, or alternatively to look up magnetic quasinormal 
modes for various magnetic held (within reasonable bounds, imposed by numerical accu¬ 
racy). As a main result, we hnd that the tensor QNMs of AdS^ Reissner-Nordstrom black 
branes can be divided into two types, referred to as “Set 1” and “Set 2”. Figure [^provides 
a schematic sketch of these two types of QNMs. Near the extremal charge, we hnd a set 
of purely imaginary modes. Set 2, which dominate the late-time behavior at larger charge 
densities. At small charge densities, these QNMs are too far down the imaginary axis to be 
seen.^ Note that also [23] show purely imaginary non-hydrodynamic QNMs in their hgure 

I. However, these modes appear in the vector huctuations (not in the tensor huctuations) 
and their behavior is not discussed in that work. Following these tensor QNMs of Set 1 
while increasing the charge density, we discover kinks in those trajectories, see e.g. hgure 
For the magnetic black brane case, we provide the values of the hrst two tensor QNMs, and 
the hrst two scalar QNMs. We hnd good agreement between our QNMs and the late-time 
behavior of the far from equilibrium system studied in [19]. These points are discussed in 
detail in Section |V| 

II. FAR FROM EQUILIBRIUM SETUPS & QUASINORMAL MODES 

The main purpose of the present paper is to provide a benchmark to compare late-time 
behavior of far from equilibrium setups to. There are in general two steps to this comparison. 
First, it needs to be determined which equilibrium state is going to be the end point of the 
equilibration process. This involves specihcation of parameters such as charge density and 

^ This fact is in stark contrast to the otherwise similar situation on AdS^ black branes where purely 
imaginary modes of electromagnetic fluctuations are present near zero frequency already for vanishing 
charge [51] . 
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magnetic field, as well as matching of the thermodynamics of the far from equilibrium system 
at late time to the thermodynamics of the equilibrium system.^ Second, at this point in 
parameter space, the solution of the far from equilibrium system has to be examined as a 
function of time.^ In general, at late times, this solution will be oscillating with a particular 
frequency and it will simultaneously decay. From this behavior the complex frequency of 
the most relevant, i.e. the lowest QNM, can be extracted.^ This kind of comparison has 
been successfully employed, for example, in |35l|36], and [T9] . 

Depending on the way in which the system is manipulated initially, it will evolve dif¬ 
ferently. This is reflected in the late-time behavior by which kinds of QNMs describe the 
approach to equilibrium. Take, for example [19], i.e. a non-equilibrium system which is 
initially sheared, between the xy-^lsne and the z-direction. This introduces a pressure 
anisotropy Ap = — (T^^) = (T^^) — (T^^) in the held theory, which corresponds holo¬ 

graphically to the metric shear huctuations {h^x — h^z) and {hyy — hzz)- Here we need to 
distinguish two cases: 

(i) The RN black brane solution enjoys invariance under SO{?)) spatial rotations when 
probed with perturbations of vanishing spatial momentum. Then {hxx — hzz) and 
{hyy — hzz) are both spin 2 tensors under these SO(3) rotations, and they obey the 
same equation of motion as the spin 2 tensor huctuation hxy. 

(ii) The magnetic brane solutions are only invariant under SO{2) spatial rotations in the 
xp-plane, because the magnetic held Fxy breaks the SO {3) symmetry. In this case 
{hxx — hzz) and {hyy — hzz) are both scalars under SO{2) rotations (even at vanishing 
momentum). 

This initial shear Ap = (T^^) — (T^^) is the kind of initial condition we have in mind 
for the main part of this paper. Other initial conditions will require diherent QNMs at late 
times. In general, one can identify the required QNMs by the symmetries which are broken 

^ For example, one needs to check if both systems are in the same phase, as characterized by the thermo¬ 
dynamic quantities. 

^ Note that the same time coordinate has to be chosen in both the far from equilibrium setup as well as 
the near-equilibrium setup for the comparison of frequencies to be meaningful. 

® Higher QNMs are also accessible after subtracting the behavior stemming from the lower ones. 
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by the initial perturbations and the background solution.® Had the initial condition broken 
a translational symmetry, for example by a gravitational shock wave initial condition [3, 
□TO], more QNMs would be excited. Then it would be interesting to consider more general 
combinations of perturbations and hence those other QNMs. We stress here, that the 
relevant helds exhibiting QNMs can be classihed according to the symmetry groups of the 
hnal equilibration state (e.g. scalars, vectors, tensors under an SO{3) rotation group). 


III. HOLOGRAPHIC SETUP 

In this section we describe a gravitational system which is dual to a particular strongly 
coupled plasma in equilibrium, namely M' = 4 Super-Yang-Mills theory at nonzero tempera¬ 
ture. Two distinct setups are discussed, one corresponding to a charged plasma, the second 
corresponding to a neutral plasma within an external magnetic held. Both of these brane 
conhgurations are solutions to the hve-dimensional Einstein-Maxwell equations. They are 
derived from the action 

S=^j d^x^[{R- 2A) - , (1) 

with K = SttGn for the hve-dimensional Newton constant Gn- A Chern-Simons term may 
be added to this action, but would inhuence neither of the black brane solutions we are 
interested in.^ 

A. Reissner-Nordstrom black branes 

Here we discuss the Anti-de Sitter Reissner-Nordstrom black brane solution to Q, which 
is dual to a charged A/" = 4 Super-Yang-Mills plasma at large N and large’t Hooft coupling 

A. 


® An additional motivation for considering the spin 2 perturbation is the fact that it is related to the shear 
viscosity through a Kubo formula |T7H55] . 

^ However, such a term potentially changes the spin 0 (scalar) and spin 1 (vector) fluctuation equations. 



1. Equilibrium solutions and thermodynamics 


The Reissner-Nordstrom black brane metric in the Poincare patch is defined by 

.2 ^2 


ds^ = — (-fdf + dx'^) + -i^dr '^, 

r^j 


with the blackening factor given by 


,, , ^ mL^ , 

/W = 1 - ^ ^ 


( 2 ) 

(3) 


This geometry has a boundary at r = oo and an outer (non-compact) horizon at some 
r = > 0. The blackening factor (|^ has six roots of which we identify the largest positive 

one with rn- At the extremal charge value q = qex two of the real roots coincide. As we 
will see later, this implies that the fluctuation equations of motion have coefficients with 
irregular singular points, while these coefficients have only regular singular points for all 
other q < q^x- The charge of the black brane couples to a f/(l)-gauge field which takes the 
form 

“ 7 ^ ’ ( 4 ) 

L 

where Q = for this to be a solution to the Einstein-Maxwell equations. 

In the rest of this paper we are going to measure quantities in units of L, which amounts 
to setting L = 1. Note that this background with m = 1 is identical to the background 
chosen in [ID] after setting the charges to zero and performing the coordinate transformation 
to the coordinate u = 

Let us collect the expressions for the relevant thermodynamic quantities in the dual field 
theory. Requiring regularity of the Euclideanized manifold at the horizon, the temperature 
and chemical potential are 


T 


2 \nrH)\ 

dvr ’ 

\/3q 


(5) 

( 6 ) 


The extremal charge is determined to be qex = derived from the condition that 

T = 0 for the extremal Reissner-Nordstrom black brane. 

It will be useful to transform and rescale our metric for the computation of the quasi¬ 
normal frequencies. We start from ([^ and perform the transformation rj^/u. In these 
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coordinates, the AdS boundary is located at m = 0, while the outer horizon is located at 
u = 1. We further rescale q —)■ qrj^. After these transformations and rescalings we get the 
metric 

ds^ = ^ {-f{u)dt^ + . , (7) 

u Au^j[u) 

with f{u) = 1 — + q'^u^, and we may use the horizon condition 0 = f{u = 1) 

mrjf = 1 + in order to remove the dependence on m. Then the metric depends only on 
two parameters namely th and q, while the blackening factor reduces to f{u) = 1 — (1 + 
q^)u^ + q^u^. ® 

We also re-write the thermodynamic quantities in these coordinates: 


T 

h 





2 - g" 1.. 



( 8 ) 

(9) 


and the extremal charge is given by q^x = Note that the relation between the original 
charge q and our rescaled q in these coordinates can be expressed in a closed form 

3/4 


q = qrjj = q 


m 


+ , 

The entropy density is given by the horizon area, and it reads 


( 10 ) 


srn — 




' H 


dvr 


dvr ’ 


( 11 ) 


where g{^) is the determinant of the spatial metric induced on the horizon. 

The energy density and pressure of the Reissner-Nordstrom black brane can be computed 
using the standard transformation to Fefferman-Graham coordinates following imiig. This 
procedure yields 


CRN 


-pRN 




( 12 ) 

(13) 


® The dependence on rn can be eliminated by a rescaling of t and x. Of course, we could alternatively 
remove all dependence on q in favor of dependence on m. 
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where k was defined in and below equation Q. Note that this implies tracelessness for the 
energy momentum tensor 




^ 3m 000^ 
1 0 m 0 0 

0 0 m 0 

0 0 0 m/ 


(14) 


such that (T^^) = 0. Even though, our ground state has nonzero temperature and chemical 
potential, both breaking the conformal symmetry, the action of A/" = 4 Super-Yang-Mills 
theory is still invariant under conformal transformations and hence operator identities are 
valid.® In particular, conformal symmetry implies tracelessness of the energy-momentum 
tensor up to contributions from the conformal anomaly given by = —aE^ — cl^ — a F‘^, 
with the anomaly coefficients a, c, a. Here, E 4 and J 4 contain various contractions of the 
Riemann curvature tensor, which vanish as our field theory lives in fiat spacetime. The third 
term —aE"^ accounts for the presence of an external field strength which also vanishes 
in the ground state discussed in this section. 


2. Fluctuations 


Having worked out the RN AdS^ black brane metric g^i,, we are now ready to introduce 
fluctuations 6 gfj,iy around this background. In this section we collect the linearized Einstein 
equations obeyed by these metric fluctuations. We choose to write these equations in mo¬ 
mentum space, after the Fourier transformation /c, m). Throughout 

this paper, we choose the radial gauge hr^ = 0 . 

Metric shear fluctuations are 2-tensors under the rotation group, and thus decouple from 
all other fluctuations. Shear fluctuations, such as h^y, satisfy an equation of motion which 
-after a field redefinition (p = = g^'^h^y- can be written as 


0 = 0 " 


f{u) -uf{u) , - f{u)k 


+ 


(15) 


uf{u) ' irjjufiuy 

We will often work with the rescaled quantities u = u/th and k = k/vn as this removes 
all factors of th from the equation of motion. This is equivalent to the coordinate rescaling 


® Because operator identities are independent of the state. 
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mentioned in footnote Note that, with the momentum in the z-direction, there exists 
one more tensor excitation given by {h^x — hyy). This combination also satishes the tensor 
equation (15) under the replacement hxy —)■ {h^x — hyy). Setting the momentum to zero 
allows two additional tensor fluctuations, namely {hxx — hzz) and {hyy — h^z)- We have 


checked that in this case all four of these tensor fluctuations satisfy equation (15) with 
k = 0, and they again decouple from all other fluctuation equations. In the present paper, 
we are interested in these sets of SO{2) tensor fluctuations hxy, {hxx — hyy) (for nonzero 
k), and S'0(3) tensor fluctuations hxy, {hxx — hyy), {hxx — h^z), {hyy — h^z) (for k = 0 ). As 
mentioned in the introduction, these tensor fluctuations correspond to pressure anisotropies 
Ap = ((T^^) — (T^^)) near equilibrium. Due to the rotational symmetry in the xp-plane, 
we know (T^^) = (T^^). At late times, a system which was initially sheared between the 
xp-plane and the 2 ; direction, should be dominated by fluctuations {hxx — hzz), {hyy — hzz) 
at vanishing momentum. For completeness we discuss here also the linearized Einstein 
equations obeyed by vector and scalar fluctuations. 

Vector fluctuations satisfy the two coupled equations 


hJ 


0 = hu" + _ 2v5«a/ , 


and 


u 


,, u {Sq'^u -2m) . 

0 — Cli H-—r-Clj + 




u 


1 ^7 


VSq [uhti + hti) 


(16) 


(17) 


/(«) • • 2f(u) 

with the gauge held huctuations a^, while indices take the values i = 1, 2, 3 at fc = 0. Using 
the constraint equation for hu, ai resulting from the rUcomponent of the Einstein equation, 
this can be transformed to a single equation 


0 = a/'+ 


+ 


1 — 2q^u 


-1 + u 1 + u 


q2y2 


„ / , - 12 g 2 (-l + u)u^{-l + u{-l + q^u)) +u\ 

' J/T, \o/T, ~oo\o ®7-ltoj 


4(—1 + u)'^u{l + u — q^v?)"^ 

Together, the vector and tensor huctuations contain all the quasinormal modes of the 
Reissner-Nordstrom black brane. Let us consider the case fc = 0 hrst. There are cou¬ 
pled equations for the scalar huctuations hu, {hxx + hyy), at- But adding and subtract¬ 
ing those equations adequately to/from each other shows that neither of these posess any 
non-trivial quasinormal mode frequencies. In other words, these scalar huctuations are 
completely determined by boundary data. In contrast to that, switching on momentum 
in the z direction, A; 7 ^ 0 , the scalar huctuations which couple to each other are given by 
htti {hxx + hyy + hzz), htz, at, a^. This set of scalar huctuations has non-trivial QNMs, for 
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example sound modes with a dispersion a; oc fc for a speed of sound Vg- We are not 
concerned with these latter QNMs in this paper and focus instead on the shear fluctuations. 


In order to find a solution to equation (15), we need to specify two boundary conditions. 


At the horizon we find the solution behaves as 


hxy = (1 — m)^ 2 { 2 - 92 ) —u) + ...] , 


(19) 


where the negative (positive) sign corresponds to the infalling (outgoing) solution. Analo¬ 
gously, we have 

Oj = (1 — 2(2-9^) I^q(o) _|_ — m) -f ...j , (20) 

for the vector modes. Quasinormal modes are dehned as those solutions to the linearized 
Einstein equations which: 

• obey the infalling boundary condition at the black brane horizon, corresponding to 


the negative sign in equation (19) and (20), and 


• satisfy a Dirichlet condition at the boundary of AdS space. 

Since these two boundary conditions are imposed at two distinct points, we need to work 
out a method to hnd solutions obeying both conditions simultaneously. This will be our 


goal in Section IV 


B. Magnetic black branes 


In contrast to the previous section, let us now consider a setup with vanishing electric 
charge, but in the presence of a constant magnetic held F^y. One such solution of action ([^ 
is the magnetic black brane [T5HT8] which is dual to an uncharged Af = 4 Super-Yang-Mills 
plasma in an external magnetic held at large N and large’t Hooft coupling A. 


1. Equilibrium solutions and thermodynamics 


For the metric and held strength we choose the following ansatz 

ds^ = -U{r)dt^+ , 

U (r) 


F = bdx A dy . 
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( 21 ) 

( 22 ) 









Then we transform to the u coordinates via r = r^u and rescale U —)■ rjjU yielding 



4:U^U{u) 


(23) 

(24) 


F = bdx A dy . 


We are interested in asymptotically AdS^ solutions to the equations of motion following 


from the action hence we set the cosmological constant A to its AdS^ value, namely 


A = —6 in units of the AdS radius L. This casts the Einstein-Maxwell equations into the 
form 


0 = 2b‘^+ 4e^^^^\u^U'{u){2V'{u)+ W'{u))+u‘^U{u){2{u{2V"{u) 

+W\u) + W\uf) + V\u){2uW\u) + 3) + ?>uV\uf) + 3W'{u)) - 3) , 

0 = 2u‘^e‘^^^^\2uU"{u) + U'{u){Au{V'{u) + W'{u)) + 3) + U{u){Au{V"{u) 

+W"{u) + W'{uf) + V'{u){4uW'{u) + 6) + AuV'iuf + 6W'{u))) - 2(6^ + , 

0 = +u\2{uU"{u) 

+U{u){4:uV"{u) + QV'{u){uV'{u) + 1))) + U'{u){8uV'{u) + 3)) - 6, 

0 = 62e-4V'H ^ 2u%U'{u){2V'{u) + W'{u)) + 2U{u)V'{u){V'{u) + 2W'{u))) - 6 . (25) 


These are the equations of motion which we solve numerically in order to determine t/, V, W 


as functions of b and of the temperature, which is associated with the horizon value of U'{u). 


All of our calculations will be performed in this coordinate system since it is convenient for 
the background. If needed, we will rescale quantities afterward in order to obtain physical 
values. 

We solve the equations of motion with the following boundary conditions at the horizon 


U = Uo + Ml(l -u) + U2{1 - uf + 0{{l - uY) , 
V = + yi(l -u) + 0{{1 - uf), 

W - U’D + 101(1 - t£) + C*((i - u)'^), 


(26) 

(27) 

(28) 


where Mq = 0 (in order to obtain a spacetime with horizon) and we can set Uq = = 0 by 
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( 29 ) 


a rescaling of the coordinates x, y, and z. The equations of motion imply 

-6 + + 9mi -6 + 56^ + 9mi 


U2 = 


Vl = 


Wi = 


12 


12 


3-62 


3mi 3mi 

]j2^-ivo 1)2 


(30) 

(31) 


6ui 6 mi 

Note that we in fact retain ten orders in the near horizon expansion in order to obtain better 
numerical solutions. However, we refrain from reproducing the lengthy expressions (for those 
M 3 , M 4 , ..., as well as M 2 , M 3 , ..., and W 2 , W 3 ,...) here. From the near-horizon expansion above 
we see that the temperature is associated with the choice of the parameter mi. The magnetic 
held will be associated with the choice of the parameter b. Choosing b and mi completely 
determines the solution to our equations of motion. 

Demanding an asymptotically AdS solution, these functions near the boundary m = 0 
behave like 


U=- 

u 


u 


B 

( 1 ) 


. 2 V 


V 

=-h 

U 


vu 


u 

B 

( 1 ) 


^ 2 W 


W 

e"" = — 
u 


H - ^ + 0 ( m °, u log u ), 

\/u 


+ 0{u^,ulogu ), 


-I- 0{u^,ulogu ), 


(32) 


where v and w are dimensionless parameters depending on the choice of Mi and b. Here, is 
the subleading coefficient in the near boundary expansion of U, which will be discussed below 


equation (39). We would like to rescale the coordinates such that the metric asymptotes to 
canonical AdS^ with a Minkowskian boundary. Hence, we rescale t ^ t = rnt, x ^ x = 
i/mx, y ^ y = y/vy, z ^ z = ^/wz, which near the boundary gives 


ds^ 


F 


—di —(^dx dy ) -|- —dz 

u u u 


+ 


dv? 


~ -dx A dy . 

V 


(33) 

(34) 


Since none of our rescalings involves the radial coordinate u, the functions U{u), V (m), W (m) 
remain unchanged. The full bulk metric after these rescalings reads 


ds‘^ = 


-Udt^+ ' ' ^'■■2^ ' ^~2 


H {dx + dy ) + 
Au^U V w 


dz , 


F = -dx A dy . 

V 


(36) 

(36) 
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In summary, in these hatted coordinates the AdS dehning function 
simply be u, giving a Minkowskian boundary metric. 


can be taken to 


The physical magnetic held, as dehned through equation (36), and the temperature in 
the boundary held theory are given by 


bA- 

V 


T = — 
271 


The entropy density is determined by the area of the horizon 


s = 


^2V+W 


47r V\/w 


1 


rn 


Attv\/w ’ 


(37) 


(38) 


where near the horizon —>■ 1 as can be seen from the expansion (26) with uq = 0 = Wq. 

In order to obtain the energy density and pressure of the dual held theory we need to work 
a little harder. We again follow nulls] in order to extract the energy momentum tensor 
from the near-boundary metric in Feherman-Graham coordinates |13|. First, we need the 
boundary expansion for the background helds in our present coordinates t, x, y, z, u: 


U = - + 
u 


u 


B 

( 1 ) 


U 


62 


u 


H -^- h + U log log u) , 


. 2 V 


-,2W 


= - + -i -^ -1- -ulogu^ + logti) , 

u 


hi) 


1(4) 
B 2 

hi) 


62 


w 


= h wu 


u 


1(4) 


Qv 


(39) 


WU 


+ 


B 2 
( 1 ) 


^^^(4) , ■M2 S/2, N 

- —u A nlogM — ^ -I- Oiu^' , logn), 

V 6 V^ 


where and are free parameters. The parameter is not free, it in fact can 


be removed by a residual diheomorphism invariance of the metric ansatz (35): sending 
M —?• m/(1 — leaves this metric invariant, and will set —)■ 0 in the expansion 

(39). In this way is seen to be related to the horizon location. 

Next, we derive the coordinate transformation which takes ns from the present coordinate 
system to Feherman-Graham coordinates which can be dehned by having a radial coordinate 
r such that Qrr = l/r2 throughout the entire spacetime, as well as the radial gauge gnx = 0 
as above. The transformation is dehned by equating the line element in Feherman-Graham 
coordinates 

ds^-PG = ^ \dr'^ + logr"^ + 0(r^))ijdz^dz^] , (40) 

to the line element in our current coordinates, i.e. 

ds^FG = , (41) 
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where ds^ is given by equation (35). Our ansatz is that the metric component coefficients 


and hO) do not depend on the Fefferman-Graham radial coordinate r, that Zq = t, zi = 
X, Z 2 = y, Z 3 = z and r is a function of only the radial u coordinate, i.e. r = r{u). Hence 
near the boundary 

r{u) = logw + 0{u^, logw), (42) 

where we will choose = 1 so that the boundary metric remains canonical Minkowskian. 


Now we plug the boundary expansions for r{u)^ equation (42), the expansions for 
U{u), Q^wiu) gjygj^ py equation (39), and the metric expansion (40), into the line 


element equation (41). Expanding the resulting equation around the boundary u = 0 and 


determining coefficients order by order we obtain the transformation for the radial coordinate 
near the boundary 


r{u) = \/u 


2 4 


u 


B 3 


U 


5/2 


log “+ O(M^ log m) . (43) 


8 "" 24r;2 

In Fefferman-Graham coordinates the near boundary metric has the expansion 


900 — —1 — ' 
gu = 922 = 1 + 

902, = l + r^ 


6 ^ — + 24:vv 


logr— + 0 {r ^), 
Zv^ 


B 

44) 


24^2 

P - 6v{vuP) + 8vuP 


logr— + 0 {r ^), 
Zv^ 


44 ) 

24'i;2 


7 2 

+ r^logr— + 0 {r^). 
Zv^ 


In this case, the holographic energy-momentum tensor is given by nam] 


(44) 

(45) 

(46) 

(47) 

where A is the renormalization scale which we will pick to be A = ^/B, and C is an ar¬ 
bitrary constant reflecting this ambiguity. We will choose C = —1/4, as this eliminates 
the P contribution to the energy density arising solely from the background magnetic held. 
The logarithmic contributions, are associated with the conformal anomaly [H]. From 




9 if - - (log(A) + C) 


equation (|47|) we obtain the energy density and pressures 

3 


e = {Too) = - 


K 


_ ,B 

■ 4 '^(4) 


Pi = P2 = (Tn) = - 


K 


4 

B 


-urA\ + — \ogB 


B^ B^ 




Pi = {Tii) = - 

K 


,B 


B^ 




(48) 

(49) 

(50) 
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FIG. 2: Energy density e of the field theory dual to the magnetic brane solution as a function 
of the dimensionless ratio (vrT)^/!?^ between the two scales in the problem, namely the boundary 
temperature T and boundary magnetic held B. At small magnetic held values, B <C the 

energy density approaches the value of the Schwarzschild AdS^ black brane, namely e = 3/4(7rT)^. 


Had we chosen a different renormalization scale A = 1,^^ this would remove all explicit 


H-dependence from the energy density, i.e. we would get (Too)a=i = —-ju^ 


k4“(4)- 


The energy density given in (48) apparently depends on both, magnetic held B and 


temperature T, seemingly independently. However, we know that the system arises from a 
conformal setting. Thus we expect that dimensionless observables should only depend on 
the ratio T'^/B of these two scales, see [19]. Therefore, it suffices to show the energy density 
(48), divided by as a function of T'^/B'^ in Fig. i Note, that our numerical data satishes 
this scaling relation fairly well, as we have checked by changing the two horizon parameters 
Ui and b. However, due to accumulating numerical errors there can be deviations of up to 
3 percent in our current data. More precisely, there is an error in our energy density that 
is on the order of 0.19 % for intermediate magnetic held (H/T^ = 0.99), and 2.3 % for large 
magnetic helds (H/T^ = 29.5). This error can be reduced by requiring higher precision. 
We refrain from pushing these deviations to zero as the current accuracy suffices for our 
purposes. The far-from equilibrium data we intend to compare to is currently known only 
up to a comparable error [T9|. 


10 


This is really A = 1/L in units of the curvature radius. 
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Note that the trace of the energy momentum tensor is given by 


(^) 


2 62 
K 2v‘^ 



( 51 ) 


This is a manifestation of the trace anomaly due to an external held F, as discussed in 
the last paragraph of Section |III A 1[ In the case at hand, that external held is simply the 


physical magnetic held B = b/v. Note that according to (51) the energy momentum tensor 
is traceless in the limit of vanishing magnetic held, as expected. 


2. Fluctuations 


In this section, we consider two types of huctuations. First, we discuss the tensor huc- 
tuations h^y and h^x — hyy, which are 2-tensors under the rotational SO{2) remaining in 
presence of the magnetic held Fxy Second, we discuss the huctuations {hxx — hzz) and 
{hyy — hzz), which are now scalars under the SO{2) and therefore couple to the other scalar 
perturbations. This set of scalar perturbations gives rise to the QNMs which describe the 
late-time behavior of the initially sheared systems we are interested in. 

Linearizing the Einstein-Maxwell equations in the huctuations, again picking radial gauge 
hry, = 0, we hnd that the mode {hxx — hyy) decouples from all others and satishes the same 
equation as the tensor mode hxy After a Fourier transformation hyy{x) oc hyy{Cj, k), 

that huctuation equation is given by 

0 = ft:. + (f - 2V + W" + 

+ ^^{U{-2b‘^e-^^ + Au^il” + 2u^ij'{AuW' + 3) - 12) 

+4M2f/2(2u(l/'2 + w" + hF'2) + 3hF') + - k‘^Ue-^^)hxy , (52) 


which, for 6 = 0, correctly reduces to the equation of a minimally coupled scalar in an AdS^ 
black brane geometry once we redehne hxy = 0/m, and choose If = l^—u, V = W = —| logw. 


We note in passing that the equation (52) at nonzero magnetic held is not simply that of a 
minimally coupled scalar held in the magnetic black brane background. 


In this section we work exclusively with hatted coordinates. So, technically, all of the coordinates t, x, y, z 
and vector/tensor components such as h^y should read t, x, y, z and respectively. In order not to 
clutter our equations, we drop the hats in this section. 
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Again, we make the ansatz h^y = (1 — u)°‘F{u), where F is regular in m at m = 1. The 
indicial exponents near the horizon are given by 


lU 


a = 


'2ui ’ 


(53) 


of which we choose the solution with the negative sign as it corresponds to infalling waves. 
We solve the fluctuation equation for h^y and require the solution to vanish at the AdS 
boundary u = 0, which gives us the tensor quasinormal frequencies for the two parameter 
family of magnetic branes (depending on the parameters b and ui). 

For k = 0, linearizing the Einstein-Maxwell equations in the fluctuations, we also hnd 
six coupled equations for the four scalar fluctuations hu, hxx, hyy, hzz- These six equations 
decouple from all other metric and gauge held huctuations, and can be found explicitly in 
the accompanying notebook | 33 ]- From them a single equation of motion for — hxx) 
can be derived. Due to the rotational symmetry in the xj/-plane we may set hyy = 0 
without loss of generality. It is convenient to work with one raised index on the metric 
perturbations, i.e. we work with hj, and h^, while hi is eliminated by making the held 


redehnition hi = x + h%. Using the tt, tu, and uu components of the Einstein equations h\ 


can be eliminated. Making use of the background equations for U, V, and W, given by (25), 
we eliminate U", U", W", and W', which also leads to cancellation of all terms involving 
hf without any radial derivative acting on it. Lastly, the xx, yy, and components of 
Einsteins can be solved for F”, hf, and a diherential equation involving only the desired 
combination of helds, x = hi — h^, which takes the form 


0 = x"(“) + a(«i‘3),\'(«) + . 


(54) 


where the coefficients a and b diverge at the horizon as usual, and depend on the radial 
coordinate u as well as on the frequency u. Their explicit form can be found in [33] . 

Finally, the scalar quasinormal modes are obtained by again demanding the Dirichlet 
condition x{u = 0) = 0 at the AdF-boundary, and the infalling boundary condition at the 
horizon 

_ 

x{u) = {1 — u) (55) 

where iL is a regular function of m at m = 1. 
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IV. COMPUTATIONAL METHODS 


In this section we review the two methods we use in order to find quasinormal modes 
(QNMs). 


A. Shooting 

The shooting method has been applied in many cases for finding quasinormal modes, 
see for example [15], HE]. Any fluctuation 0 (e.g. a metric component, gauge field, scalar 
fluctuation) generically satisfies an equation of motion which is of second order in the radial 


coordinate, see equations (15) and (52) for our specific cases at hand. Hence, one has to 
specify two boundary conditions in order to numerically solve these equations. The basic 
idea of the shooting method is to specify these two boundary conditions at the horizon, 
and then adjust the free parameter given by the frequency u to find the desired solution 
at the boundary. One of these two horizon boundary conditions is the infalling boundary 


condition. This amounts to picking the negative sign in equations (19) and (53). The 


remaining boundary condition is a mere normalization coefficient. In order to ensure that 
the solution is a quasinormal mode one has to vary the frequency u until the field vanishes 
= 0 at the boundary r = tb- This can be efficiently achieved using a numerical 
optimization procedure such as Mathematica’s FindRoot [ 32 ], nesting this FindRoot with 
the NDSolve that finds the numerical solution, see [33] for more details. It is now clear 
that the shooting method depends crucially on the convergence properties of the horizon 
expansion and the precision of the initial values provided from evaluating that horizon 
expansion numerically at a numerical cut-off r = rcut-off near the true horizon r^- 

Once a QNM is found at a particular point in parameter space, e.g. specified by q, we 
can easily obtain QNMs for an e-neighborhood around this initial point in parameter space, 
i.e. q + e. We merely give the previous QNM frequency to FindRoot as a best estimate of 
the result, and find the corresponding QNM at that neighboring point in parameter space. 
Generally, the performance increases when e is decreased. However, there is only a limited 
region in parameter space which is accessible to this method without increasing the working 

Fluctuations can couple to each other, see [niQE] for a systematic method for finding quasinormal modes 
of such coupled systems of fluctuation equations. 
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precision, or improving performance of the algorithm otherwise, as we discuss in the next 
paragraph. 

One advantage of the shooting method is that it can be used for hnding quasinormal 
modes of spacetimes that are only known numerically. While the method is also quite fast, 
it still has its limitations. One such limitation is that performance decreases rapidly as 
quasinormal frequencies with large imaginary part (Imw ^ T) are considered. See |15] for 
a discussion and example calculations for this and related issues. In order to overcome this 
limitation one has to compute more coefficients of the near horizon expansion in order to 
be able to provide initial values further away from the horizon location r^j (without leaving 
the radius of convergence of the horizon expansion). 


B. Continued fractions 


The continued fraction method has also been put to great use in numerically determining 
quasinormal modes, as in [I1I19]. This method makes use of an elegant mathematical theorem 
originally due to Pincherle, and later generalized in [50]. Instead of dealing directly with 


the differential equation of motion, i.e. (15), this technique works with the power series 


solution expanded about a singular point, such as the horizon. Generically, the radius of 
convergence of such an expansion is limited by the distance to the next nearest singular point 
of the differential equation. On the other hand, Pincherle’s theorem, and its generalization, 
give a criterion to determine when this radius of convergence is increased. 

An increased radius of convergence allows the calculation of quasinormal modes for the 


following reason. The fluctuation’s equation of motion (15) has a regular singular point at 


the boundary, where it’s characteristic exponents are A± = ± 2 . This means that a solution 
there will generically behave as 


hxy = + aiT + ■ ■ ■ + a 4 r ^ log r + ■ ■ ■ 


+ bor + bir H-, 


-3 


(56) 


where Oq and bo are the two free constants determining the expansion and, importantly, 
since the characteristic exponents differ by an integer, the logarithmic term starting with 
04 is needed so that the two solutions are linearly independent. The other Oj and 6 * are 
determined in terms of oq and bo by the equation of motion, and crucially 04 vanishes only 
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when ao is zero. Therefore, by applying Pincherle’s theorem to the power series of the 
infalling solution at the horizon, we can determine a criterion on the frequency ca such that 
this series is convergent at the boundaryConverging to an analytic function, it cannot 
have the logarithmic term 04 , which in turn means that this solution has oq = 0. Such 
u therefore determine which infalling solutions have the leading behavior hxy ~ as 

r —)■ 00 , that is, they are Dirichlet at the boundary: they are the quasinormal frequencies. 

1. Pincherle’s theorem 

Pincherle’s theorem, and its generalization, relates a recursion relation, such as that 
coming from a power series solution of a differential equation about a singular point, to a 
continued fraction. The convergence of the continued fraction corresponds to the existence 
of a minimal solution of the recursion relation, where we recall that a sequence hn is minimal 
if 

lim — = 0 , 

for all other sequences Qn satisfying the recursion relation. A minimal solution implies an 
increased radius of convergence of the power series solution, beyond the next singular point. 
For a full proof see the reference [50] , the important results are collected in Appendix 

The continued fraction method of calculating the quasinormal frequencies is useful as a 
numerical check of other methods. Beyond that it is powerful because it is often faster than 
the shooting method, and can explore further into the complex plane with greater accuracy. 
Its major shortcoming is that it can only be applied to hnd the quasinormal modes of 
fluctuations around analytically know backgrounds, and for that reason is inapplicable to 
the magnetic branes also of concern at present. 

2. Application to Reissner-Nordstrom 

In applying the continued fraction method to the calculation of quasinormal modes, 
the hrst thing needed is an expansion of the equation of motion about a singular point. 

This method fails when the boundary is not the next nearest singular point to the horizon. 
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The desired singular point is the horizon, as that allows the infalling nature of the per¬ 
turbation to be made directly in the ansatz. We will start with the background metric 
of the Reissner-Nordstrom black brane in the coordinates of Q. The perturbation ansatz 
Qxy = uhxy{u) is made and the Einstein-Maxwell equations are linearized in hxy, as 

It is useful to work with the rescaled quantities Cj = oj/th and k = k/rn 


in Section 


III A 2 


as this eliminates th from the equation of motion. One last change of radial coordinate is 
made: u = 1 — r(2 — This is done as the horizon is now at fjj = 0, which allows a 
cleaner series expansion there. The charge dependent rescaling simply helps with numerical 
stability in the region of extremality, q —)■ y/2. 

The equation of motion for h^y has regular singular points at the horizon and the bound¬ 
ary = 1/(2—as well as at the radial coordinates^^ r± = {—l+2q‘^±^yl + 4q‘^)/{2q‘^{2— 
g^)^). Determining the leading behavior near the horizon, the characteristic exponents are 
a = ±tCj/{2{2 — g^)), with the lower sign corresponding to the desired infalling condition. 
Peeling off this factor, i.e. hxy{r) = r°f{r), the equation of motion can be written as 


0 = 


/('■) 'll 


(57) 


i=0 


i=0 


i=l 


where r*, s,, and ti are coeffcients that depend on g, u, and k and can be found in the 
notebook 1^. 


Postulating a series solution, /(r) = ^ Qr*, the coefficients q must obey the recursion 

i=0 

relation 

5 

ri i + {n- i)Si + {n - i){n - i - 


Cr?, — 


-E 

2=1 


Cn—i 


with the initial conditions 


Co — 1, Cl — 


CoCo 


C2 = 


nso + n{n — l)fi 


ci(ro -h si) -h Con 


(58) 


C3 = 


C4 = 


Sq 2so + 2 ti 

C 2 (ro -F 2si -F 2^2) + Ci(ri -F S2) + Cor 2 
3so + Qti 

C 3 (co + 3 si -|- 6^2) + C2(ri -|- 2s2 + 2^3) -|- Ci(r2 -f- S3) -|- Cor3 

4so + 12ti 


(59) 


The location r_|_ is always further from the horizon than the boundary is, but at g > •\/3/2 the singular 
point r_ is the next nearest singular point to the horizon. This corresponds to charges greater than 
91.7%qex for the unsealed charge q, and invalidates the continued fraction method in this region, as 
mentioned above. 
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Comparing to Parusnikov’s Theorem 2, see Appendix this 6 term recursion relation 
implies m = 4, and the 5-dimensional vector pn needed to construct the continued fraction 
has components 


Pi,n 


+ {n + i - 6)s6-i + {n + i - Q){n + i - 7)t7_i 
nso + n{n — l)ti 


( 60 ) 


The desired continued fraction /o is now dehned by the limit given in (A2). In practice, 
one calculates for large f; changing this parameter allows one to understand the conver¬ 
gence properties of the limit, and obtain any desired accuracy. The hnal step in calculating 
the quasinormal frequencies is to demand that a minimal solution to the recursion relation 
(58) exists by hnding the values of u such that Theorem 7/8 of Appendix [A| is obeyed for 


the initial conditions (59). Explicitly, we hnd the solutions of 


C4 = 





( 61 ) 


which is a high order polynomial in 0. Note that /s is dehned in the appendix (below 
Theorem 7/8). Figurej^shows the movement of the hrst three quasinormal modes calculated 
in this way as the number of iterations i is increased from 3 to 41, in steps of two. The modes 
generically move up and to the right, and converge rather quickly. In practice, 50 iterations 
were used to obtain the 10“® accuracy in the continued fraction data quoted throughout. 


V. RESULTS 

In this section we describe the quasinormal modes (QNMs) which are present in the 
tensor and vector huctuations around Reissner-Nordstrom black branes, and in the ten¬ 
sor huctuations around magnetic black branes. In the latter case we also study the scalar 
quasinormal modes. We discuss the behavior of the quasinormal modes as a function of 
the relevant parameters and provide explicit data for various points in parameter space for 
comparison. Throughout the present section we employ different units: we report quasinor¬ 
mal frequencies in units of mass density, by which we mean that we display values in tables 
and hgures which are computed at the parameter value m = 1, if not specihed otherwise. 
Where indicated, we display the frequencies in units of temperature (ttT), or of mass den¬ 
sity after a rescaling with the inverse horizon radius, u ^ uj/th (stemming from a rescaling 
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FIG. 3: Convergence of the continued fraction method. The first three quasinormal modes (circles, 
squares, diamonds) at g = /c = 0 are plotted in the complex frequency plane as the number of 
iterations is increased. 

of the time coordinate t —?■ rjjt in the metric ([^).^^ The latter values are referred to as 
“rescaled QNMs”. The charge density is given either in units of mass density, i.e. q in units 
of (m/3)^/^, or in dimensionless values denoted by q. Momenta k are given in units of mass 
density m, momenta k and k are the rescaled versions. 

A. Quasinormal frequencies of Reissner-Nordstrom black branes 

We hnd two distinct sets of modes, see the schematic plot in hgure The hrst set (set 
1) also exists at zero charge density and has been analyzed in pQ. The second set of modes 
(set 2) consists of purely imaginary quasinormal modes. We observe the latter to approach 
the origin at and above values of g = O.dqex- These purely imaginary modes of set 2 are not 
visible at all at small or zero charge densities, as they are lying far down the imaginary axis. 

In agreement with previous results [1], we hnd a tower of quasinormal modes with fre¬ 
quencies that have nonzero real and imaginary part. The values obtained from our continued 
fraction method are in very good agreement with those obtained from the shooting method 

Note that in the coordinates discussed in this sentence the energy density is now temperature and charge 
dependent. 
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mode # 

shooting 

cont. fraction and Starinets 

1 

±3.1194516 -i2.7466757 

±3.119452 -12.746676 

2 

±5.1695210 -i4.7635701 

±5.169521 -14.763570 

3 

±7.1879308 -i6.7695650 

±7.187931 -16.769565 

4 

±9.1971992 -^8.7724814 

±9.197199-18.772481 

5 

±11.2026887 - H0.7740258 

±11.202676-110.774162 


TABLE I: Consistency with previous data. The values of the hve lowest lying quasinormal mode 
frequencies u for the tensor fluctuations, e.g. hxy, are shown in units of mass density. The first 
column labels the mode number n = 1, 2, 3, 4, 5. The second column shows the quasinormal mode 
frequencies obtained from the shooting method, the third column shows those obtained from our 
continued fraction method which is in complete agreement with the values previously obtained by 
Starinets in [T]. These results are computed at vanishing charge density. 

(8 significant digits agree), and both are in good agreement with previons results found by 
Starinets. These results are computed at vanishing charge density. This agreement is a non¬ 
trivial result, as our continued fraction calculation is quite different from that of Starinets [T]: 
ours is a 4 dimensional vectorial continued fraction, while [1] uses a (1 dimensional) scalar 
continued fraction. See table HI 

At increasing charge density and hxed momentum k = 0, the modes of set 1 move as 
indicated in £gure|^ displayed in the frequency u. There is an obvious kink in the trajectories 
at larger charge densities q ~ O.Gqex- However, plotting the QNMs in the rescaled frequency 
variable a), no kinks are visible. Only upon multiplying a) with the horizon radius, those 
kinks appear in the QNM trajectories. These kinks appear suspiciously close to the region 
where convergence of the horizon expansion within our shooting method gets increasingly 
worse.However, we have checked that the location of the kinks does not depend on the 
numerical horizon and boundary cut-offs that we choose. Independently, the continued 
fraction method conhrms the kinks in the trajectories at the same locations. That said, 
when plotting the frequencies in units of temperature (vrT), or in units of we observe no 

This is caused by fact that the inner horizon, located at r_, is approaching the outer horzion, located at 
ru- The kink in the quasinormal mode trajectories appears where the difference — r_| is comparable 
to the horizon cut-off located at Vcut-off^ i-e. {th — r_| « |r// — Vcut-offl- 
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q/qex 

Cj (rescaled) 

U}/{7tT) 

0 

±3.11945 - i2.74668 

±3.11945 -i2.74668 

0.1 

±3.11652 - i2.75222 

±3.12256 -i2.75756 

0.2 

±3.10759-^2.76953 

±3.13227-^2.79153 

0.3 

±3.09236 - i2.80090 

±3.14994-i2.85305 

0.4 

±3.07057-^2.85111 

±3.17859 -^2.95141 

0.5 

±3.04312 - i2.92958 

±3.22533 - i3.10499 

0.6 

±3.01749 - i3.05496 

±3.31032 -^3.35142 

0.7 

±3.03304 - i3.25409 

±3.50619 - i3.76173 

0.8 

±3.16481 - i3.47084 

±3.9918 - i4.3778 

0.9 

±3.32234 - i3.68357 

±5.0231 - i5.56924 


TABLE II: Reissner-Nordstrom QNMs of hxy at q ^ 0, k = 0. Lowest QNM frequency in set 1 
for increasing charge densities q (in units of (m/3)^/^) given in fractions of the extremal charge 
density qex- The left column provides the QNM frequencies rescaled with a factor rj{, while the 
right column shows the frequencies of the same modes in units of (vrT). 


kinks whatsoever, merely smooth trajectories. Hence, these kinks should be regarded as an 
artifact of expressing frequencies in units of mass density. This point is explored in more 
depth in one of the accompanying notebooks [5H] . 

At vanishing momentum, k = 0, the leading QNM of the metric shear fluctuation hxy 
is listed for increasing charge q of the Reissner-Nordstrom black brane in table [TTj The 
right column shows the frequency values in units of (vrT), while the left column shows 
the rescaled values which we obtain directly from our numerical procedure because therein 
we have worked with conveniently rescaled quantities, as seen from the metric ([^. Those 
rescaled coordinates are convenient for our near-equilibrium setup. For a comparison to the 
late-time behavior of physical quantities obtained from numerical non-equilibrium codes the 
values in the right column are more suitable. In comparison with ra. our lowest QNM 
values agree up to deviations of 0.2 percent, as can be seen from Table 1 in [T9] . 

Fixing the momentum at larger values k = 1/2 and k = 1 does not change the qualitative 


behavior of the modes, as seen from tables III and IV, respectively. Also in the quasinormal 


mode trajectories computed at hxed nonzero momentum k = 1/2 we observe kinks, as 
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FIG. 4: Reissner-Nordstrom QNMs of h^y at q 0, k = 0. The three trajectories of the three 
lowest quasinormal modes are shown as the momentum is fixed to k = 0 and the charge is increased 
from q = 0 to q = 0.9 qex in increments of Aq = Qex/dO. The uppermost point in each trajectory 
(top right point) corresponds to vanishing charge q = 0. For increasing charge, the modes each 
move to more negative imaginary frequency values. 

before at A; = 0. See figure However, at large fixed momentum, k = 5, we observe 
that the kinks disappear from the regime accessible to our numerical methods, as seen 
in hgure Another interesting exercise is to £x the charge density and then follow the 
quasinormal modes while increasing the momentum. The result of this exercise can be seen 
in hgure for hxed charges q = 0, O.bqex, 0.8qex (blue, red, brown curves, respectively). At 
large momentum, the trajectories all asymptote to the same “attractor” curve, regardless of 
the charge. This seems reasonable because at large momentum the scale set by the charge 
is negligible compared to the scale set by the momentum, in other words all trajectories at 
large momentum should asymptote to the zero charge trajectory, which is conhrmed by our 
data. 

Now we turn to the second set of quasinormal modes (set 2), namely the purely imaginary 
ones mentioned before.The member of this set 2, which is closest to the origin a; = 0, is 

We are grateful to Michal Heller for pointing us to these kinds of modes. Further, we are grateful to 

Julian Sonner for discussions on the interpretation of these modes. 
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mode ff 

1 

2 

3 

q = 0 

±3.1779741 -i2.7290839 

±5.2112434-14.7519066 

±7.2214147-16.7604134 

Q — 0.25 Qex 

±3.1413918-z2.7474627 

±5.1413864-14.7936400 

±7.1185609-16.8281811 

q = 0.5 qex 

±3.0262710 - z2.8281461 

±4.9367320 - 14.9919741 

±6.8493201 -17.1646578 

q = 0.75 qex 

±2.9252877 - z3.0924526 

±4.9163769 - 15.3955189 

±6.8500839 - 17.6457490 


TABLE III: Nonzero momentum. Reissner-Nordstrom black brane quasinormal mode frequencies 
u for the tensor fluctuations, e.g. hxy, at fixed momentum k = lj2 and increasing charge density 

q- 


mode ff 

1 

2 

3 

q = 0 

±3.3440382 - 12.6808260 

±5.3317107-14.7178383 

±7.3188762 - 16.7331804 

Q — 0.25 Qex 

±3.3102870 - 12.6965339 

±5.2651313 - 14.7562293 

±7.2195106-16.7971758 

q = 0.5 qex 

±3.2045995 - 12.7642384 

±5.0692976 - 14.9339576 

±6.9544811 -17.1066419 

q = 0.75 qex 

±3.0864530 - 12.9801352 

±5.0079788 - 15.3207002 

±6.9204927 - 17.5887749 


TABLE IV: Nonzero momentum. Reissner-Nordstrom black brane quasinormal mode frequencies 
Lo for the tensor fluctuations, e.g. h^y, at fixed momentum k = 1 and increasing charge density q. 
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EIG. 5: Reissner-Nordstrom QNMs of hxy at q ^ 0, k = 1. The trajectories of the three lowest 
quasinormal modes are shown as the momentum is fixed to A; = 1 and the charge is increased from 
q = 0 to q = 0.9 qex in increments of Aq = qex/^0. 
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FIG. 6: Reissner-Nordstrom QNMs of hxy at q ^ 0, k = 5. The trajectories of the three lowest 
quasinormal modes are shown as the momentum is fixed to A: = 5 and the charge is increased from 
q = 0 to q = 0.9 qex in increments of Aq = gei/40. 
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FIG. 7: The trajectory of the three lowest quasinormal modes are shown as the charge is fixed to 
q = 0,0.bqex, or O.Sqex while the momentum is increased from A: = 0 to A: = 20 in increments of 
Ak = 0.1. 

shown along with the lowest QNM from set 1 (and its complex conjugate) in hgure That 
hgure shows the trajectories of those three poles as the charge is increased from q = 0.6 
to 1.0 in steps of 1/40. Note that we have chosen to display the rescaled frequencies here 
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(stemming from a rescaling of the time coordinate in the metric ([^, as discussed before) d® 
While the lowest mode known from set 1 moves (slowly) away from the real axis, the lowest 
mode of set 2 move on the negative imaginary axis towards the real axis. At a value of 
^crossing ~ 0.655 the purely imaginary mode and the lowest mode of set 1 have identical 
imaginary parts. See hgure At charge densities q > qcrossing ~ 0.655, the lowest purely 
imaginary mode is the lowest of all QNMs. 

As the charge density approaches its extremal value from below, the modes of set 2 are 
pushed towards the origin and closer towards each other. We conjecture that this may be 
understood as the onset of a branch cut forming along the negative imaginary frequency 
axis, in analogy to the near-extremal AdS^ case |5lll52] Our calculations indicate that 
the modes of set 1 do not approach the imaginary axis in the extremal limit. All of the 
QNMs of set 1 and set 2 remain in the lower half of the complex frequency plane. In other 
words, the fluctuation h^y does not seem to cause any instabilities in the extremal limit (this 
statement is limited by numerical accuracy). 

We have confirmed selected QNMs of set 2 with the continued fraction method to an 
accuracy of 4 digits.Various checks of both of our numerical methods indicate that these 
purely imaginary modes of set 2 are a robust feature which is not an artifact of the numerical 
methods. Consequently, the QNMs of set 2 are dominating the late-time physics of the 
dual held theory at fairly sizable temperatures, namely at charges q > ^crossing ~ 0.655, 
corresponding to temperatures T < Tcrossing ~ 0.228. 

Our lowest tensor QNM frequency from set 1 agrees very well with the QNMs extracted 
from the late-time behavior of the fully time dependent setup at the corresponding charge 
densities. Table 1 in [19] shows a maximum deviation of 0.2% at fairly large charge densities. 
At charge densities much smaller than the extremal one, we hnd much better agreement. 
Note that Table 1 in [19] only probes charge densities for which the lowest QNM of set 1 

We also find it convenient to express the charge density in its dimensionless form here, i.e. we are using 
q as opposed to q. 

Just like in that case, at extremality, i.e. in the zero temperature limit, conformal symmetry is still broken 
by the nonzero chemical potential. 

High accuracy at all parameter values is difficult to achieve in the purely imaginary QNMs. The problem 
is that, at particular values of the charge g, the zeros of (j){u = 0, w; q) move through poles of the function 
(j){u = 0,Cj;q). Those poles are located at w = in{q^ — 2) (with integer n), where the continued fraction 
method is known to yield false QNMs [U 02] . 
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FIG. 8: At large charge densities, a purely imaginary QNM enters the low frequency regime. 
We display the rescaled frequencies after oj —)• oj/rn- The trajectories of the three quasinormal 
modes closest to the origin a; = 0 are shown as the charge increases from q = 0.6 to 1.0 in 
increments of Ag = 1/40. The momentum is fixed to A; = 0. The black dots mark the points of 
the respective trajectories with q = 0.6. The lowest QNM in set 1 and its complex conjugate mode 
move downward as the charge increases (blue downward-pointing triangles). The lowest QNM in 
set 2 moves upward as the charge increases (red upward-pointing triangles). 


still dominates the late-time behavior. The authors of [19] use the charge density q which 
takes its extremal value at ~ 0.6204 with the mass parameter m = 1, while 

the extremal charge value in the tilded coordinates is g = a/2. The relation between the 
two charges is g = g(m/(l -|- g^))^'^^. Therefore, 90% of extremal charge in [19] is identical 
to 58% of extremal charge in the tilded coordinates. The authors of [19] report numerical 
difficulties with larger charge density values, whereas it is straightforward for us to reach 
more than 70% of extremal charge in the tilded coordinates. 

For comparison to [23], here we also compute the vector modes. The two lowest vector 
QNMs (members of Set 1) at vanishing momentum k = 0 are given for various charges 
in table |V} Figure IT shows their trajectories with increasing charge density, and we have 
chosen the same units for our frequency as the authors of 123, i.e. we have plotted u/ (27rT). 
Note that the authors of [25] display results at nonzero momentum only, while all of our 
vector QNM results are at vanishing momentum. Comparing to hgure 1 from [25], we 
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FIG. 9: The imaginary part of the lowest QNM in set 1 is shown (line with negative slope) as 
a function of charge q. This is compared to the imaginary part of the purely imaginary QNM 
which is the lowest member of set 2 (line with positive slope). The two lines intersect at the charge 
^crossing ~ 0.655 with an imaginary part of Imw ss —3.50. 

observe qualitative agreement to the extent allowed by the visual comparison We compare 
to the left plot in hgure 1 of [2S]- light of our own results, we interpret their green 
data set (r_/r_|_ = 0.5) as similar to our vector QNMs at intermediate charge values, i.e. 
roughly half of the extremal value of q. In that case, there is already one mode of Set 2 
visible, sitting on the imaginary axis in their hgure 1. We interpret their blue data set 
(r_/r+ = 0.95) as showing only modes of Set 2 whereas the modes of Set 1 are outside 
the plot range of their hgure 1. Our numerical data at vanishing momentum does not 
show the hydrodynamic (dihusion) mode closest to uj/{2tiT) = 0 indicated by the box in 
hgure 1 of [25]. At vanishing charge g = 0, we recover the values, u:/{2tiT) = ?7,(±1 — i) 
with n = 0, 1,..., already obtained analytically in [53]. Just like for the tensor huctuation, 
also in this vector huctuation a second set of purely imaginary QNM (Set 2) wanders up 
the imaginary axis. And like the purely imaginary tensor QNMs, at a critical charge value 
also these imaginary vector QNMs become more relevant to the late time behavior than the 

The figure in [25] shows vector QNMs at nonzero spatial momentum. An analysis of the vector modes 
at nonzero momentum is beyond the scope of this work, hence we content ourselves with a qualitative 
comparison. 
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mode ^ 

q/qex = 0 

q/Q ex — 0*8 

q/qex = 0.9 

Set 1 (complex QNMs) 

1 

±1.0000000-il.OOOOOOO 

±3.3687082 -il.5042378 

±6.3036022 - i2.9733044 

2 

±2.0000000 - i2.0000000 

±5.8977388 - i5.4108751 

±9.0680850 -i6.1472415 

Set 2 (purely imaginary QNMs) 

1* 

not accessible 

±0 - i2.5337209 

±0 - i2.3875204 

2* 

not accessible 

±0 - i3.8897030 

±0 - i3.4700708 


TABLE V: RN AdS^ vector QNMs. The values of the two lowest lying complex and two lowest 
lying purely imaginary quasinormal mode frequencies in units of temperature, i.e. uj/{2ttT), for 
the vector fluctuations are shown at selected charge densities q. 

modes of Set 1. In summary, the vector channel QNMs show exactly the same structure 
and qualitative behavior as the tensor QNMs. 


B. Quasinormal frequencies of magnetic black branes 


Here we report quasinormal mode frequencies of the magnetic black branes discussed in 


Section III B On the field theory side this setup allows us to examine the late-time behavior 


of neutral strongly coupled plasma in an external magnetic held being sheared and then 
relaxing. The shear occurs between the a:?/-plane and the z-direction. On the gravity side 
we have introduced scalar huctuations hu^ hxx, hyy, h^z at vanishing momentum, which have 


common QNMs determined by the huctuation equation (54). We also study the QNMs of 


the tensor huctuation hxy as discussed in Section IIIB2 obeying equation (52). 


Let us begin with the scalar QNMs, see Fig. 11 We provide a collection of sample 
values for the two lowest QNMs at three distinct points in (T, i3)-parameter space, see 


table |VI[ The lowest of the scalar QNMs in Tab. |V^ is in reasonable agreement with the 
QNM values extracted from the late-time behavior in |T9], which are compared to our 


results in Tab. VII The hrst column shows seven distinct values of B/T^] we conhrmed 


that the (dimensionless) value of the energy density over agrees at these points within 
a few percent deviation between the non-equilibrium setup [12] and our magnetic brane 


thermodynamics from Sec. Ill B Agreement between the thermodynamic values for the two 
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FIG. 10: The two lowest complex and two lowest purely imaginary RN AdS^ vector QNM 
frequencies in units of temperature, i.e. uj/{2ttT), are shown for increasing charge densities q/qex = 
0.800,0.805,0.810,..., 0.895, 0.900. Like the tensor QNMs, also these vector QNMs fall into two 
sets. Set 1 (blue down triangles) moves down in the complex plane as the charge density is increased. 
Set 2 (red up triangles) moves up along the imaginary axis as charge density is increased. The 
black dots mark the QNMs at q/qex = 0.800. 


mode # 

R/r2 = 0 

R/r2 = 12.953 

R/r2 = 30.161 

1 

±3.1194506-^2.7466751 

±3.4079677 - i 3.0777933 

±3.8221728-i 3.5553514 

2 

±5.1700747-i 4.7637826 

±5.4444069-i 5.5805876 

±5.6661971 - i 6.6938458 


TABLE VI: Magnetic brane scalar QNMs. The values of the two lowest lying quasinormal mode 
frequencies u for the scalar fluctuations, e.g. hxx — hzz, are shown in units of mass density. The 
columns show the quasinormal mode frequencies for three distinct values of B/T^. In the B —)■ 0- 
limit, these two modes connect smoothly to the two lowest modes of set 1 of the RN black brane 
at vanishing charge density g = 0, i.e. to the ones found earlier by Starinets pQ. 


setups is best at small values of B/T'^ (about 0.2 percent) at point 2, then deviations rise 
for bigger values, and asume their maximal value (about 2.3 percent) at point 7, i.e. for 
relatively large B/T"^ = 30.161. 

Now we turn to the tensor QNMs at vanishing momentum, two of which are shown in 


Fig. 12 Note that the lowest tensor QNM approaches the imaginary frequency axis at 
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parameter point 

QNM 

Cj 

late-time [TH] 

UJ 

relative deviation [in %] 

ARew, Almti; 

1 (B/r2 = 0.000) 

±3.119450641 - i 2.746675059 

±3.1195-i 2.7466 

±0.00002, -0.00003 

2 (B/r^ = 0.990) 

±3.122662648 - i 2.750477496 

±3.124-i 2.73 

±0.04283, -0.74451 

3 (B/r2 = 5.344) 

±3.197822654 - i 2.838566372 

±3.217-i 2.79 

±0.59970, -1.71095 

4 (B/r2 = 12.953) 

±3.407967707 - i 3.077793253 

±3.480-f 2.96 

±2.11364, -3.8272 

5 (B/T2 = 17.821) 

±3.538016365 - i 3.224040043 

±3.634-i 3.09 

±2.71292, -4.15752 

6 (B/T2 = 22.836) 

±3.661199784 - i 3.364602583 

±3.780-i 3.21 

±3.24484, -4.59497 

7 (B/T2 = 30.161) 

±3.822172759 - i 3.555351392 

±3.98-i3.38 

±4.12925,-4.93204 


TABLE VII: Magnetic brane scalar QNMs versus far-from equilibrium late-time oscillations. The 
value of the lowest lying quasinormal mode frequency u for the scalar fluctuations, e.g. hxx — hzz, 
are compared to the oscillation frequency extracted from the late-time behavior of the full time 
evolution in m- 



Re 

FIG. 11: Trajectories of the two lowest QNMs of scalar (hxx — hzz) fluctuations around magnetic 
black branes: The magnetic field is increased from B/T‘^ = 0 (at an energy density of e/T^ = 73.4) 
to 15.5 (at an energy density of e/T^ = 147.6). These values are generated by fixing the near 
horizon parameter ui = 2, and choosing 6 = 0 to 19/20 in steps of Ab = 1/20. Both, real and 
imaginary part increase when the magnetic field (and with it the energy density) is increased. The 
second QNM frequency moves a much greater distance in the complex frequency plane than the 
first. 
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mode jf 

B/T2 = 

5.342 

B/T2 = 22.80 

B/r2 = 30.00 

1 

±3.0104279 - 

i2.8338442 

±2.1829303 -i3.4535515 

±1.8555798, -i3.7287196 

2 

±4.9749665 - 

f4.8923780 

±3.6619677 - f5.3641251 

±3.6224530 - i5.5936925 


TABLE VIII: Magnetic brane tensor QNMs. The values of the two lowest lying quasinormal mode 
frequencies uj for the tensor fluctuations, e.g. h^y, are shown in units of mass density. The columns 
show the quasinormal mode frequencies for three distinct points in (T, B) parameter space. In the 
B —7- 0-limit, these two modes connect smoothly to the two lowest modes of set 1 of the RN black 
brane at vanishing charge density g = 0, i.e. to the ones found earlier by Starinets [I]. 



FIG. 12: QNMs of hxy fluctuations around magnetic black branes. The two lowest QNMs (both 
members of set 1 in the limit of vanishing magnetic field, see Fig. at vanishing momentum are 
shown as the magnetic field parameter b is changed from 6 = 0 to 6 = 50/40 in steps of Ab = 1/40. 
These modes both move down and to the left with increasing b. This covers a magnetic field 
over temperature range of BjT"^ = 0.0001, ...,34.4555, which corresponds to energy densities of 
e/T^ = 73.4, ...,397.4. We have fixed the near horizon temperature parameter to ui = 2. 

large magnetic fields B S> 10. For the second QNM, however, our numerical results are 
inconclusive at those large magnetic field values. We provide a collection of sample values 
for the two lowest QNMs (from set 1) at three distinct points in (T, i3)-parameter space in 
Tab.IVnH 
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VI. CONCLUSION 


Strongly coupled plasmas can be studied as they thermalize, because we have various 
holographic methods at our disposal. What is the correct picture for such an evolution from 
a ferocious state of matter to a tame thermal equilibrium? At late times in the thermalization 
process, we at least have a fairly clear picture on the gravity side: just like the overtones 
of a bell that was gently struck and then left to ring down freely, various quasinormal 
modes (QNMs) of the black brane are excited. All but the lowest lying QNMs (with the 
smallest imaginary part) have died down at late times because all but the lowest ones are 
strongly damped. This picture is valid for small deviations from equilibrium and at late 
times. However, one may be inclined to think of intermediate states, and possibly even the 
initial state of the gravitational setup, as a superposition of (not necessarily quasinormal) 
modes.Then, it would be interesting to know which of the modes are excited at what 
time during the evolution, what energy they carry, and how quickly they dissipate. Is all 
the energy hrst transferred from the low frequency modes into the high frequency modes, or 
vice versa? Our analysis in comparison with the far from equilibrium data from [12] allows 
a few speculations on which we comment below. 

This present paper demonstrates the computation of tensor fluctuation {hxy) quasinor¬ 
mal modes within two distinct black brane systems asymptoting to AdS^: First, Reissner- 
Nordstrom (charged) black branes. Second, magnetic black branes (uncharged). In the 
latter, we also compute the QNMs of scalar fluctuations h^, hyy, h^z, while in the 
Reissner-Nordstrom case we re-compute the vector QNMs originally found in [25]. We pro¬ 
vide sample values in tables ill |III[ |IV| , [y] |VI[ |VII[VIII[ Two notebooks are provided 
online [33] for computation of QNMs at user-selected charge density q and momentum kz, 
or for reference QNM frequency values at particular magnetic helds B, and temperatures T. 
All these quasinormal mode frequencies should be used in order to compare to holographic 


A detailed analysis |54j found that a quasinormal mode decomposition approximates the full nonlinear 
time evolution surprisingly well for all tested initial conditions. In order to support the idea of formulating 
time-dependent systems in terms of a well-defined set of modes, we refer to the analysis of AdS^ black hole 
QNMs performed in global AdS, and then reinterpreted in the Poincare patch [^. There, the authors 
relate that QNM spectrum to a time-dependent held theory stressing similarities to heavy ion collisions. 
Note also the perturbative weak held calculation of a gravitational collapse geometry in presence of 
charge [56]. 
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calculations of anisotropic fluids far from equilibrium. The late-time behavior of the latter 
should be governed by oscillations and decay that is well approximated by particular QNMs. 

As one example, we have compared our QNMs to a system [12] which starts with a metric 
that is anisotropic between the x|/-plane, and the z-direction, and which is then observed as 
it freely evolves either in the presence of (i) charges or (ii) a magnetic held on the black brane. 
We hnd good agreement. In the case of a charged plasma, dual to the Reissner-Nordstrom 
AdS^ black brane (case (i)), our lowest tensor mode QNM (associated with huctuations h^y 
or equivalently h^x — h^z) agree with the values provided in [12] up to 0.2% deviation at 
large charge densities.In the case of a magnetic plasma, dual to the magnetic AdS^ black 
brane case (case (ii)), our lowest scalar QNM (associated with huctuations hxx — hyy) is in 


good (at worst 5% deviation) agreement with [T2|, see Tab. VL These hndings conhrm our 
reasoning from Section [IT] that the late-time behavior of a relaxing system is determined by 
the symmetries of the background it relaxes towards, and by the symmetries of the initial 
conditions. In the Reissner-Nordstrom case, the initial shear in the setup [12] (far from 
equilibrium) is a tensor under SO{2) rotations in the x|/-plane. At late times, only the 
corresponding tensor QNM matters. In the magnetic black brane case, on the other hand, 
the same initial shear is a scalar. At late times, only the corresponding scalar QNMs matter. 

Furthermore, we provide QNM values for Reissner-Nordstrom black branes at nonzero 
momentum. We may think of these QNMs as being useful for comparison to far from 
equilibrium systems which started out with an initial condition depending on the z-direction. 
Colliding shock waves with a 2 ;-dependent prohle may serve as an example for this kind of 
system. 

At sizable charge densities, a novel set of QNMs dominates the late-time behavior of 
the charged black brane system, and hence also the late-time behavior of the dual charged 
plasma. In addition to the well known set of QNMs (set 1) of the Reissner-Nordstrom black 
brane, we also found a set of purely imaginary QNMs (set 2) as schematically shown in 
Fig.ll These modes enter the low energy spectrum at and above a particular charge density 
q > 0.655, corresponding to a sizable temperature T < Tcrossing ~ 0.228, with the lowest 
member of that set becoming the new lowest lying QNM of the sytem. As the sytem is 


This may have to do with the lowest mode of set 2, beginning to play a significant role at these larger 
densities; but not yet dominating the late-time behavior. 
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pushed toward extremality, based on our numerics, we conjecture that these modes of set 2 
are going to be pushed towards each other and form a branch cut. An analogous behavior 
is seen in the AdS^ case mm- 

In our second system, the magnetic black branes, we also hnd tensor QNMs at vanishing 


momentum. We provide sample values in table |VIII| , and provide a notebook with samples 
of QNMs at various temperatures, and magnetic helds. 

It is interesting to note that all of the QNMs of magnetic branes move deeper into the 
complex frequency plane as the magnetic held is increased. This should be interpreted as 
a faster dissipation of each mode after it has been excited. It is then tempting to argue 
for a faster equilibration of the non-equilibrium system if all the QNMs decay faster at 
larger magnetic helds. However, this does not seem to be the case |I9]. Instead, (charge 
or) magnetic held individually seem to have a rather negligible ehect on the equilibration 
time [19]. This may be interpreted in (at least) two ways: as evidence for the fact that 


(i) the standard QNMs (which we computed here) are not a good measure for equilibration 
of this system while it is far from equilibrium, and one needs to hnd other ways of 
characterizing this process, or 


(ii) only the lowest of the QNMs plays a signihcant role in the equilibration process. This 
is a possibility, because the lowest magnetic brane QNM remains virtually static, 
while the higher magnetic brane QNMs change signihcantly with the magnetic held. 


as illustrated in Fig. 11 


We point out a fundamentally diherent behavior for a subset of the electric black brane 
QNMs. As noted in the previous paragraph for magnetic brane QNMs, most of the electric 
brane QNMs move deeper into the complex plane at increasing charge density. The notewor¬ 
thy exception are the QNMs of set 2, see Fig. [Tj Those QNMs move closer to the origin of 
the complex frequency plane while the charge is increased (and simultaneously temperature 
approaches zero). This behavior indicates a decreased dissipation of these modes, i.e. 
these modes become more and more long-lived. It is tempting to ask if these modes become 


Numerically, we can only verify this up to a certain lowest charge value, at which these imaginary QNMs 
are simply too deep in the complex frequency for our methods to find them at the current accuracy. 
However, all available results indicate that these imaginary QNMs move closer to cc = 0 with increasing 
charge, at any given charge value. 
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part of an effective field theory description near zero temperature. Such a description could 
contain these long-lived modes whose frequency is much smaller than the charge density, 
i.e. uj ^ q, near zero temperature. This would be in anology to the standard hydrodynamic 
description of long-lived modes which satisfy the relation uj T. Such an effective held 
theory description may be useful for the study of black holes near extremality. It should be 
stressed that these modes are not propagating. 

To provide a brief outlook, it would be interesting to consider the other huctuations of our 
two systems. All huctuations could be studied at nonzero momentum in the ^-direction. The 
same can be done for a nonzero momentum in the ^/-direction. In that case, one would then 
compute the QNMs of metric huctuations which are vectors or scalars under the remaining 
rotation group. Among those QNMs one expects dihusion and sound modes, i.e. modes 
with dispersion relation a; oc /c^ and u o: k, respectively. This would then allow one to 
compare to the late-time behavior of more complicated far from equilibrium systems. 

Finally, we note that the shooting and continued fraction methods take quite a long 
time to hnd the present results. More realistic, and thereby more complicated systems, do 
encourage the development of improved methods for the computation of QNMs. 
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Appendix A: Generalization of Pincherle’s theorem and continued fraction method 

In its generalization, Pincherle’s theorem relates m + 2 term recursion relations to m- 
dimensional continued fractions. Such vectorial continued fractions require use of a pro¬ 
jective space, all needed notions will be dehned as they arise. The hrst useful result is 
Parusnikov’s Theorem 2 which provides the correspondence 
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Theorem 2: 

i) interuption-free m-dimensional continued fraction; 

fo = bo+ g , 

b2 + -- 

m+1 

iii) m + 2 term recursion relation: Qn = Pm+ 2 -j,ngn-j, 

i=i 

where: Pn = (l/on, ' ,bm,n/cin)', geuerically Vij denotes the i-th component of 

vector Vj; and “interuption-free” means that an is non-zero. 

The construction of the m-dimensional continued fraction is at hrst opaque, as it requires 
“dividing by a vector,” but it can be simply accomplished by the following procedure. First, 
given the vector p„, dehne 


U ■■■ 0 
1 ■■■ 0 



Pl,n 

P2,n 


Qn — Po ■ ■ ■ Pn, J 


yO • • • 1 Pm+l,n J 


/o 1 0 ■■■ o\ 

0 0 1 ■■■ 0 

0 0 0 ■■■ 1 

\i 0 0 ■■■ 0/ 


(Al) 


These matrices act on m -|- 1-dimensional projective vectors, denoted (/i : ■ ■ ■ : fm+i)- 
For the map [(/i : ■ ■ ■ : fm+i)] = (/i//m+i, • • ■ ,/m//m+i), the m-dimensional continued 
fraction can be computed 


fo = lim = lim a ^3 ^Q^e^+i, (A2) 

where e^+i = (0 : ■ ■ ■ : 0 : 1). Theorem 2 therefore provides a way to construct an m- 
dimensional continued fraction given an m -|- 2 term recursion relation. 

The next theorem provides a powerful criterion for when the continued fraction /q con¬ 
verges. 


Theorem 6: 

An m-dimensional continued fraction fg converges. 


There exists a minimal sequence hn that satisfies the recursion relation given in Theorem 2. 


Parusnikov’s Theorems 7 and 8 now give a recursion relation that allows an explicit 
determination of hn- 
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Theorem 7/8: 

The minimal basis obeys: 


hn — 


-1 


fm, 


n+1 


m—1 


hn—m ^ ^ 


-k 


/>+i 




m,n+l 


hn—m+j 1 


(A3) 


where /„ is the continued fraction fn = hn + a„/(6„+i + ■■■). This formula allows the 
calculation of the quasinormal frequencies, as explicitly demonstrated in section IVB2[ By 
setting n = m, the hi become the hrst m+1 coefficients of the recursion relation generated by 


the equation of motion. One then uses (A2) to plug in the components of fm+i as a function 
of a). Demanding the veracity of the formula in Theorem 7/8 determines the frequencies that 
give the minimal solution, and hence are the quasinormal frequencies, as discussed above in 
section IIV B 11 
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